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Artificial Modification of the Earth’s 
Radiation Belt*‘ 


S. F. Singert 


; Summary 


The fundamental theory of the behavior of charged 
particles, such as electrons and protons, in a magnetic 
field was set out by the Norwegian astrophysicist 
Stormer, and then developed more fully by the Swedish 
astrophysicist Alfvén. His theoretical ideas were ap- 
plied by us to a hypothesis that charged particles from 
the sun and from cosmic rays could be trapped in the 

.earth’s magnetic field. This hypothesis was developed 
originally to explain the origin of magnetic storms and 
aurora. It received striking support by Van Allen’s 
- discovery of the “radiation belt’? which must consist 
of such trapped particles. From theory we can calculate 
the length of time for which the particles stay trapped 
and even make suggestions as to their origin, but to 
‘get greater precision in theoretical predictions it is 
_ essential to do experiments. One of the most interesting 
experiments which suggests itself is that of injecting 
artificially accelerated particles directly into the radia- 
tion belt region. One promising method would consist 
of carrying a small electron accelerator aloft in either 
a satellite or even a high altitude rocket. The most 
direct place for conducting the experiment is near the 
equator at an altitude of from 600 to 1000 miles. The 
particles which are injected, preferably electrons of 
the order of 2 Mev, would rapidly diffuse in a type of 
“napkin ring” region around the earth. This diffusion 
only requires a few seconds since the particles move 
with nearly the speed of light. Electrons then stay 
trapped in this napkin ring region, diffusing out of it 
very slowly. The primary reason for their outward 
diffusion may be collisions with the very rarified atmos- 
phere which is still present even at these high altitudes. 
Therefore by studying the number of remaining elec- 
trons from day to day we can learn more about the 
density of the atmosphere in a region of space where 
densities can normally not be measured. The simplest 
way to determine this density is by firing sounding 
rockets daily through this napkin ring region and 
determining the electron intensity profile from day to 
day and see how it decays. We would expect such a 
region to last for a few days to a few weeks, depending 
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on the altitude and depending on what actual air densi- 
ties exist at these altitudes. 

The injection apparatus would only weigh of the 
order of 500 pounds including its power supply, and 
could provide us with a very substantial increase in the 
radiation belt intensity in a small region. The reason 
for this striking increase with such small equipment 
is again due to the fact that the earth’s magnetic field 
has such good trapping properties. The sounding rockets 
which are sent up after the injection to check up on the 
fate of the electrons can be quite simple and need only 
to carry Geiger counters with associated circuits weigh- 
ing on the order of a pound or less. A project of this 
type is therefore feasible immediately. 


Historical Introduction 


The earth’s magnetic field may to a good approxi- 
mation be represented by that of a dipole. More than 
fifty years ago Carl Stérmer, a Norwegian astro- 
physicist, became interested in the problem of the 
motion of a charged particle in such a dipole field 
(1). He developeda very elaborate theory which may be 
summarized as follows: A charged particle starting off 
far from the earth may or may not reach a point under 
consideration. In general, there will be some directions 
of incidence which are allowed. These directions are 
therefore contained within the ‘“‘allowed cone” while 
the rest of the solid angle may be called the ‘‘forbidden 
cone”’ since particles of a certain energy coming from 
infinity cannot reach the point under consideration in 
that direction. (Of course, for particles of very low 
energy all directions are ‘‘forbidden,’”’ while for particles 
of extremely high energy, which travel through the 
field without much deflection, all directions are 
‘“allowed.’’) 

Stérmer undertook this study in order to explain 
the origin of the aurora, assuming that they were 
produced by particles originating at the sun and then 
entering the earth’s magnetic field. However, it turned 
out that his studies were more appropriate to the be- 
havior of the much higher energy cosmic rays. St6rmer 
and the workers following him mostly paid attention 
to those particles which could reach the earth from 
infinity. But what if we trace back the trajectory of a 
particle along a forbidden direction? It turns out that 
these particles could not have come from infinity, but 
that their orbits are “trapped,” i.e., confined to the 
vicinity of the earth. Stérmer, in fact, showed that 
there are regions around the earth in which particles 
could move but that these regions have no physical 
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connection to infinity, and he concluded therefore that 
they are of no interest (Figure 1). 
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Fie. 1. Diagram showing the allowed and forbidden spaces 
in the vicinity of the Earth’s dipole according to Stormer. 
Small changes of the dipole field can open up to infinity the 
allowed region in the immediate vicinity of the dipole; this 
region extends narrow funnels into the Earth’s atmosphere in 
the auroral zones; for particles of very low energies, encoun- 
tered in solar corpuscular streams, the simpler peturbation 
methods of Alfven can be used to discuss the problem of 
trapping in the Earth’s dipole field. 


First to consider the possibility that there might 
be particles in trapped orbits seems to have been 
Hannes Alfvén, a Swedish astrophysicist, well known 
for his development of the science of magnetohydro- 
dynamics. Alfvén investigated the orbits of cosmic rays 
in the vicinity of the sun, assuming the sun to have a 
dipole field. He showed by simple considerations that 
this dipole field could not be perfect because the earth 
possessed a magnetic field which would perturb the 
sun’s field near the earth. Therefore the earth could 
scatter particles into these forbidden regions and the 
particles might remain there for long periods of time 
(2). This idea was developed in greater detail and with 
great elegance by Professor John Wheeler of Princeton 
University and his students (3). However, the idea 
was not taken up in a practical way because further 
work showed that the sun probably did not have a 
dipole field and that in any case the motion of the 
conducting interplanetary gas would disturb the sun’s 
magnetic field considerably. 

In the meantime, laboratory experiments had been 
performed which clearly showed the existence of par- 
ticles in the forbidden regions. Birkeland, in Norway, 
(4) at the end of the last century was first to develop 
the potentialities of the “terrella,” a model earth having 
a magnetic dipole field. This was followed many years 
later by very precise work in Alfvén’s Institute by 
Malmfors (5), Brunberg and Dattner (6), and by 
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Block (7), leading to model experiments on the trajec- 
tories of cosmic rays and on the origin of the aurora. 
Willard Bennett (8) at the Naval Research Laboratory 
carried out an especially elegant demonstration of the 
trajectories of electron streams in the vicinity of a 
terrella. His photographs of the ‘“Stérmertron” show 
very clearly the existence of particles in forbidden 
regions around a dipole. 

The existence of trapped particles in the earth’s 
magnetic field (9) seems to have been considered first 
in 1956. Even though a single particle could not enter 
into the forbidden region (as shown by Stérmer’s 
theory), Singer reasoned that if a large number of 


particles were to arrive from the sun, their collective — 


action could perturb the strict dipole field sufficiently 
to allow entry into the trapping regions. 


By perturbing the magnetic field, particles are 


scattered into the normally forbidden regions so that 
after a solar corpuscular cloud has passed the earth, a 
certain number of particles may remain in this region. 
The motion of the particles was discussed and it was 
shown that they would not immediately hit the earth 
but, instead, remain trapped for long periods of time 
(9) and would be removed only by interactions with 
the very tenuous atmosphere or by the scattering effects 
of magnetic field inhomogeneities. It was also suggested 


that cosmic ray albedo could make a contribution to — 


this trapped radiation (9). 

It must be understood that the assumption of par- 
ticles in the forbidden regions was purely a hypothesis 
and was developed mainly to explain magnetic storms 
and aurora. Subsequently, Singer developed a more 
detailed theory for the main phase of magnetic storms 
which makes use of the drift of these trapped particles 
to explain the formation of a ring current (10). But 
since the occurrence of a ring current was not certain 
and since in any case its mechanism was not under- 
stood, the hypothesis of the existence of trapped par- 
ticles in the earth’s dipole field was still unsupported. 

An important advance from an experimental point 
of view was made by Van Allen and his colleagues who 
measured for the first time the incidence of charged 
particles in the auroral zone (11). Those particles were 
presumed to be electrons. Although the existence of 
particles had been suspected for a long time as a source 
of the aurora, this constituted the first direct measure- 
ment. Now if the hypothesis of trapped particles was 
correct, then one should be able to see these particles 
also at lower latitudes but at higher altitudes, namely 
along the same line of force (see Figure 2). Accordingly, 
a proposal was made to detect auroral particles by 
means of a high altitude rocket launched at low lati- 
tudes (12). This experiment was one of the main ob- 
jectives of the U.S. Air Force FAR SIDE project (13); 
an extreme high-altitude probe which was launched in 
October 1957 from Eniwetok and penetrated well into 
the radiation belt region. 

Van Allen’s dramatic discovery of a layer or belt 


Fig. 2. Schematic indication of the reflection and trapping 
of a charged particle by the earth’s magnetic field. The par- 
ticle’s velocity makes a pitch angle a» to the magnetic line of 
force in the equatorial plane. As the particle moves into a 
region of stronger magnetic field (lines closer together), its 
pitch angle increases until it reaches 90° at which point the 
particle is reflected and moves into the opposite hemisphere 
where it will again be reflected from a similar turning point. 
With a smaller value of a» the particle will approach closer 
to the earth’s surface. For extremely small values of a the 
particle can enter the earth’s atmosphere. With ay = 90° 
the particle remains in the equatorial plane. Both soft and 
hard radiation belt particles are trapped in this manner. 


of radiation surrounding the earth (14) (which became 
evident through a Geiger counter in the first Explorer 
satellite) gives us now direct evidence for the existence 
of trapped particles in the vicinity of the earth. 

In what follows we will be discussing a summary of 
the trapped particle theory and some of its applica- 
tions, in particular the possibility of modifying the 
radiation belt. We want to discuss in detail the possi- 
bility of injecting particles into the earth’s magnetic 
field in order to study the storage properties of the 
magnetic field and in order to measure the atmospheric 
density at very high altitudes. 


Theoretical Treatment: The 


Particles 


The theoretical study of the radiation belt revolves 
on the question of how particles may enter or be released 
into the trapping regions and in what manner and for 
how long the particles move inside this region. 

Turning to the second problem, the motion of par- 


Trapping of 


ticles in the magnetic field, we find here interesting 
developments reaching back to nearly 100 years. A 
magnetic dipole field, such as that of the earth, can be 
represented as follows. The equations below give the 
magnetic field strength B as a function of distance r 
from the dipole center and as a function of the polar 
angle } measured from the equatorial plane. 


= = (1 + 3 sin’ \)"? (1) 


where M = 8.1 X 10” gauss-cm’, the earth’s dipole 
moment. 

The lines of force which indicate everywhere the 
direction of the field are then given by an equation 
of form 


r = 7. cos (2) 


where r, is the distance from the dipole center where 
the line of force intersects the equatorial plane. Pro- 
ceeding along a line of force starting from the equatorial 
plane it can be seen by inspection that the lines of force 
crowd together as one approaches the earth and that 
therefore the field becomes increasingly stronger. It is 
not possible to express the magnetic field strength B as a 
function of distance along the line of force in closed 
form; however, various approximate formulas can be 
developed to describe B (15). In polar coordinates, 
along a given line of force 


_M(1+3sin'))” 


B 
3 cos? 


(3) 


A charged particle moving in a magnetic field B with 
velocity V experiences the well-known Lorentz force 
(e/c)(V X B). If we decompose its velocity vector V 
into a velocity Vy parallel to B and V\ perpendicular 
to B we will find the particle moving along a helical 
trajectory whose pitch angle a is given’ by 
tana = V,/V,. The particle thus advances along the 
field with velocity V,, meanwhile gyrating about the 
line of force. The period and radius of gyration are 
given by 


27mc 

te me eB ) (4) 
mV Cc 

pais (5) 


It was first realized by the French mathematician 
Poincaré (16) that a charged particle which moves in 
an inhomogeneous magnetic field and approaches a 
region of increasing strength of the magnetic field can 
be reflected back into the region of weaker field, i.e., 
that V,) decreases to zero and changes its sign. Stérmer 
(1) carried out many numerical calculations of particle 
orbits which clearly exhibit this feature of reflection by 
the stronger magnetic field. However, it was Alfvén 
who first showed that this property of the motion of a 
particle could be understood in terms of an approxi- 
mate conservation law (17), namely the conservation 
of the particle’s magnetic moment yu, defined as the 
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current generated by its gyration times the area en- 
closed. Thus 


B= 


a a > Si 
- a = smV mV? sin a (6) 


a1 


Alfvén showed that this quantity remains constant to a 
good approximation in a magnetic field which varies 
smoothly (17). If the magnetic field possessed in- 
tensity variations which are severe over a distance 
shorter than the particle’s radius of curvature, then the 
magnetic moment of the particle would no longer stay 
constant; i.e., when 


psrad By/.B > -1. Ci 


This latter case of course holds for the motion of cosmic 
rays in the earth’s dipole field since their radius of 
curvature can be very large. As a matter of fact, no 
proper theory has as yet been developed to indicate 
precisely how the magnetic moment may vary in a 
dipole field and how rapidly this approximate con- 
stancy is lost. Numerical calculations indicate that the 
magnetic moment is conserved even upon many re- 
flections to a surprising degree (18), but for long life 
times such questions have to be considered with much 
more precision. However, in what follows we have as- 
sumed strict conservation of the magnetic moment. 

Alfvén also pointed out (17) that time variations of 
the magnetic field with a time scale shorter than the 
gyration time of the particle could also change its 
magnetic moment; 1.e., for 

te a8h (55 > he (8) 
dt 
We shall not consider these effects in the development 
of the theory but return to them only later. 

Keeping in mind that the kinetic energy of a particle 
must remain constant in its motion in a time-invariant 
magnetic field we see that sin’ a/B must remain con- 
stant. It can therefore be seen that sin @ reaches its 
maximum value of 1 (and a becomes 90° and Vy be- 
comes zero) at a certain value of the magnetic field B. 
The particle cannot move beyond this point since sin a 
cannot increase beyond 1; therefore the particle will be 
reflected at this turning point. The location of the 
turning point depends in the first instance on the par- 
ticular line of force along which the particle is moving 
and secondly on its initial pitch angle ao (defined as 
the pitch angle of the particle when it crosses the 
equatorial plane of dipole). 

We have therefore a picture of particles spiraling 
along lines of force from the northern into the southern 
hemisphere to be reflected there, returning to the 
northern hemisphere, being reflected at a symmetrical 
turning point, and therefore remaining trapped in the 
earth’s magnetic field. Particles with very small pitch 
angle can of course reach much lower altitudes and may 
therefore escape into the dense atmosphere almost im- 
mediately (see Fig. 2). 
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Superimposed on this oscillating motion is a slo 
azimuthal drift of the particles which may give rise to 
current. The drift velocity is given by 


Va = kV p/r (9 ) 


where k is of the order of 1, and depends on the valu 
of ao. In fact we were led to the hypothesis of trapped 
particles in order to explain through this drift the ring: 
current responsible for magnetic storms-(10). As of this” 
date the existence of such a current has not been verified” 
experimentally although the satellite observations cer- : 
tainly verify the existence of trapped particles, and the 
drift. associated with them follows from the theory. 


Lifetimes of Particles in Trapped Orbits : 

In order to find the lifetime of particles in these 
trapped orbits (and from this their intensity), we must 
discuss how particles once trapped can be released. The 
physical ideas were again discussed in earlier publica-_ 
tions (10) and are basically connected with a change in 
pitch angle. If the pitch angle of a particle can be made 
smaller then it will travel further along a line of force, 
and reach deeper into the dense atmosphere; having 
gone deeply into the atmosphere, it will lose enough 
energy so that it will no longer be effective. The pitch 
angle can be changed basically by two methods (10).. 
First, by collisions with atoms or ions of the outermost 
atmosphere of the earth, and secondly, by scattering 
from discontinuities and short-scale time-variations of 
the earth’s magnetic field (as indicated by equations 
(7) and (8)). Which of these two processes is most 
important depends on the particular situation which 
we discuss. For example, for low altitudes at the equator, 
the first process is likely to be most important since at 
low altitudes the atmospheric densities are still quite 
high. Also near the equator the earth’s magnetic field 
is “quiet”? as compared to the state of the field in the 
auroral region. In the auroral region at high latitudes, 
and along the same line of force at high altitudes, the 
magnetic field is more lilely to be disturbed by the con- 
tinuing incidence of solar corpuscular streams and by 
the turbulence contained in these streams. In fact, 
bunching of the trapped particles may themselves 
originate magnetic field fluctuations which propagate as 
hydromagnetic waves and serve to accelerate some of 
the particles while releasing them through the change 
of their pitch angle. While this phenomenon may be of 
importance to the explanation of the aurora (19) At 
will not concern us in this paper. 

Before discussing any theory for the origin of the 
radiation belt, it.is important to calculate the lifetime 
of particles in trapped orbits. This can be done with 
some precision for particles at low altitudes near the 
equator. There, as indicated earlier, the residual at- 
mosphere is likely to provide the most important 
mechanism for limiting the particle’s lifetime. We have 
investigated this point very carefully and found that 
relativistic electrons (i.e., electrons with energy greater 


than the electron rest mass of 0.5 Mev) are scattered 
‘in interactions with atoms, ions and other electrons, 
and that this angular scattering produces a type of 
“wandering” in pitch angle. As shown in Ref, (20), 
this “random walk” in pitch angle can be described by 
a diffusion equation which leads then to a means for 
estimating the rate of escape of the particles from the 
trap. On the other hand, for protons having an energy 
of about 100 Mev the main process is found to be 
energy loss by collisions (21). It has been shown that 
this will decrease their energy to a point where the 
protons are no longer detectable and therefore are no 
longer of experimental interest. 


Origin, Intensity, and Altitude Distribution of 
the Radiation Belt 


_ When the high intensity radiation was discovered in 
Explorer I satellite, it was at first thought that this 
radiation might be due to the trapped particles re- 
“sponsible for the magnetic storms and aurorae. At high 
latitudes or large altitudes solar corpuscular radiation 
should certainly be the predominant source. However, 
it is difficult to see how this radiation can diffuse to 
low altitudes at the equator; at least no quantitative 
theory has as yet been given for such diffusion. There- 
fore, we have investigated other injection mechanisms. 
We have worked out in some detail a possible mecha- 
nism based on cosmic ray albedo neutrons, which could 
provide trapped particles in the equatorial and low 
latitude regions (21). The same idea has been sug- 
gested independently in the Soviet Union by Lebedin- 
sky and Vernov (22); and by Rothwell, Gold (23) and 
Kellogg (24). Kellogg’s calculations, however, give 
values for lifetimes and intensities which differ con- 
siderably from our values. 

The albedo neutron injection works as follows: 
Primary cosmic rays smash into the earth’s atmosphere 
in the vicinity of 20 to 25 km (70,000 feet) and dis- 
integrate atmospheric nuclei. Neutrons are released 
with energies up to several 100 Mev; many of them 
travel upwards and are therefore called “cosmic ray 
albedo.” 

Being uncharged, the neutrons travel along straight 
lines right through the earth’s magnetic field out into 
space, but free neutrons are radioactive. A small 
fraction of the fast neutrons, about 10 per cm or 
about one in a million in 100 miles, change into protons 
and electrons while traveling through the field. Com- 
pared to the protons which have an energy of the order 
of the neutron energy, the lifetime of the electrons is 
rather small. The protons immediately spiral around 
the line of force, a fraction having a small pitch angle 
are lost immediately as they get too deep into the at- 
mosphere, but the rest are trapped. 

The essential part of the calculation can be under- 
stood as follows: The concentration N of trapped par- 
ticles follows from a balance between the injection rate 
Q and the rate of removal; the latter can be expressed 


in terms of a lifetime 7 and is found to be determined 


mainly by collisions with atmospheric atoms and ions. 
Roughly then, we find 


N per cm’ ~ Q (per em® — sec) X T' (sec) (10) 


But the experimentally determined quantity is the flux 
J, given by J (perem” — sec) ~ N (per em*) X 
V (cm/sec). A typical injection rate is 10 */em*-sec 
(for protons near 50 Mev) ie., 1 particle is injected 
into a given em® of space every 300,000 years. But the 
theoretically calculated “lifetime” is so long that the 
particle concentration builds up to about 10‘ per em® 
at 1000 km (600 miles).* With a velocity of ~10” 
cm/sec the flux is ~1000. 
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Fie. 8. Radiation belt intensity and atmospheric density 
and composition at the equator. The atmospheric density was 
derived from satellite data, the low altitude values from the 
atmospheric drag, and a high altitude point at 1900 km from 
the radiation belt intensity and the theory of Ref. 21. The 
exosphere critical level is found at 530 km; above this level 
collisions between atoms can be neglected and the temperature 
remains constant at 1500° K. Oxygen, because of its smaller 
scale height falls off very rapidly with altitude until at about 
1000 km hydrogen predominates. This region then we call 
the “‘outer exosphere”’. Using these new values for the atmos- 
pheric density the radiation belt intensity can be calculated 
as a function of altitude and is shown in the graph. 


The lifetime which is calculated is found to be in- 
versely related to the atmospheric density. Since the 
latter decreases with altitude, the lifetime and there- 
fore intensity of the radiation increase with altitude. 
The radiation flux (in particles per cm’-second-stera- 
dian) is found to vary with atmospheric density d as 
J ~ 10 “/d at low altitudes near the equator (21). 
With a radiation belt flux of 3000 at 1900 km given by 
Explorer IV (22) we can therefore derive the density of 
the atmosphere at 1900 km as 2.5 X 10° gram per 
em’, i.e., an H-atom concentration of about 1.5 X 10° 


* The situation is analogous to a water tank which has a 
very tiny leak; then a small injection rate can build up a high 
level. 
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per cm®. Combining with density values at lower alti- 
tudes from satellite orbit data we can construct a model 
of the exosphere. Its temperature is 1500° K. The criti- 
cal level is found to be at 530 kilometers and transition 
to the ‘outer exosphere” (where neutral hydrogen pre- 
dominates) occurs at 1000 kilometers. (See Figure 3). 

Using these values for the atmospheric density calcu- 
lated from the one observed radiation belt point, we 
can now calculate also the variation of radiation in- 
tensity with altitude. For the cosmic ray albedo radi- 
ation belt a maximum is reached at the equator at 
about 14 earth radii (near 6000 miles) beyond which 
its intensity drops. The radiation belt which is of solar 
origin and connected with aurora should reach its 
maximum near 5-6 earth radii (20,000—24,000 miles). 
The belt’s intensity ultimately decreases because the 
trapping properties of the field deteriorate, and dis- 
appears when the interplanetary magnetic field takes 
over, probably near 10 earth radii. These theoretical 
predictions can be checked experimentally; the moon 
shots should give an excellent opportunity for such 
measurements. 


Experimental Results 


Experimental results on the radiation belt have been 
obtained by Van Allen and his colleagues in Explorer 
satellites (14) and by S. N. Vernov, A. E. Chudakov 
and their colleagues in Sputnik IIT (26). (In a different 
Sputnik III experiment, V. I. Krasovsky (27) has 
measured large fluxes of electrons and photons, es- 
pecially in the auroral zone, and has pointed out the 
great importance of this radiation to the heating of the 
upper atmosphere in the auroral zone.) The experi- 
mental results have as yet not been presented in a 
definitive way. 

As of this date the nature of the radiation belt par- 
ticles has not actually been determined. However, in 
our view, there are in reality two radiation belts co- 
existing at the same time: (1) the hard radiation belt, 
just discussed, which is of cosmic ray origin, and can be 
calculated with great accuracy from cosmic ray data 
and theory. (2) In addition we also have the soft 
radiation belt which is of solar origin and connected 
with magnetic storms and aurorae. However, its lo- 
cation is different; it surrounds the hard radiation belt 
like a halo, although there must be a region of overlap 
(see Figure 4). As a consequence the radiation belt 
picture changes at different latitudes and different 
altitudes. Table 1 indicates schematically the proper- 
ties of the hard and soft radiation belt of the earth. 

The energy spectrum to be anticipated is shown in 
Figure 5. The high energy portion is the hard radiation 
belt, the low energy portion is due to the soft radiation 
belt. The intensity of the soft belt will vary depending 
on location and time; for example, it will be absent at 
low altitudes near the equator and will be a maximum 
in the auroral zone following intense solar activity. As 
a consequence of this complicated mixture of spectra, 
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Fic. 4. Distribution of radiation belts around the earth, 
derived from the theory in Ref. 21 and 30. The lines emanating 
from the earth indicate the lines of force of the earth’s mag- 
netic field. The closed lines give contours of equal radiation 
intensity for the hard radiation belt, which is of cosmic ray 
origin. Its maximum is in the vicinity of 1.5 earth radii. Its 
intensity begins to increase at about 600 km (according to 
the theory) and may reach out as far as 10 earth radii (40,000 
miles) at the equator. Note that the regions around the poles 
are clear. The soft (auroral) radiation belt is indicated by 
cross hatching and overlaps somewhat with the hard radiation 
belt at higher latitudes and larger altitudes. The maximum of 
the auroral particle radiation belt should occur between 4 to 6 
earth radii and depending on the detector and shielding used 
may exceed in intensity the hard radiation belt. 


experiments with thin absorbers are likely to give in- 
consistent results. Since the hard radiation belt should 
be quite constant with time, the use of shields having a 
thickness of about 10 g-em~ or more, should give re- 
producible flux values. 


Artificial Modification of the Radiation Belt 


One may take advantage of the excellent trapping 
properties of the earth’s field in two different ways. If 
the lifetime is as long as one calculates based on at- 
mospheric scattering only, then it may be possible to 
substantially decrease the intensity level in the belt by 
providing additional absorbing bodies such as large 
satellites (28). 

However, before carrying out any large-scale experi- 
ments, one would like to study the trapping properties 
of the field experimentally. For this purpose, it would 
be most interesting to inject particles into the magnetic 
field directly at high altitudes and trap them (28). 
Electrons are better suited for the injection than 
protons since they are much easier to generate and since 
their radius of curvature is very small. For this reason, 
they become quite insensitive to magnetic field fluctu- 
ations and to questions of the constancy of the magnetic 
moment. (Compare equations (7) and (8)). For ex- 
ample, a 2 Mev electron has a radius of curvature at an 
altitude of 1000 km above the equator, which is only 


3 km as compared to 70 km for a 100 Mev proton. The- 


injection is best carried out by means of a linear ac- 
celerator which is carried aloft in a satellite ; even a 
high altitude rocket may be adequate. 

It turns out that an injection time of only a few 


TABLE 1 


Properties of Earth’s Radiation Belts 


Hard (cosmic-ray)* Soft (Auroral) 

- Origin From the decay of cosmic | From Solar Corpuscu- 
ray albedo neutrons lar streams with 
which come out of the subsequent _accel- 
earth’s atmosphere. eration near the 

earth. 
Nature Protons between 10-400 | Electrons and_ pro- 
Mev tons of less than 1 
Mev 

Shielding | Difficult to absorb and | Easily absorbed 

shield. 
Location | Equator and low latitudes | Auroral latitudes 
Close to the earth Far from the earth 
Time de- | Constant in time Variable, increases 
pend- when sun is active 
ence 


' ® From Reference (21). 
b From Reference (10). 


minutes is quite sufficient to give an increase of the 
electron intensity observable over the natural radiation 
belt intensity (28). 
- The way in which such an experiment might work is 
as follows: The accelerator should be as small and as 
light as possible. This suggests a linear accelerator of 
the type that uses a klystron as an electron source and 
also as a means of supplying energy to a wave guide 
which accelerates the electrons. The instrument can 
be constructed to be extremely light-weight since 
the vacuum system can of course be dispensed with. 
This leaves only the electron gun, the tube structure, 
focusing magnet, and wave guide. By supplying a small 
magnet the electron beam can be doubled up and the 
whole accelerator can be folded into a smaller space. 
The power supply requirements can be estimated 
reasonably as about 50 kW for a few minutes. Applying 
a normal weight factor for batteries, we derive a total 
payload weight of about 500 pounds which is sur- 
prisingly small compared to what one might expect 
from laboratory types of electron accelerators. 
The mean output current of such an accelerator 
would be about 500 microamperes giving therefore 
3 X 10” electrons per second. If we inject over an 
altitude interval of 10 km, then the volume into which 
the electrons are released and diffuse into is of the order 
of 10° cm*. Its shape resembles a “napkin ring.’’* The 
mean injection rate thus becomes 3 X 10° electrons per 
em’ per second. For simplicity we may assume that the 
injection is isotropic, i.e., the satellite or rocket is 
tumbling so that electrons are injected into all di- 
rections (it should be noticed here that it is not neces- 
sary to provide a well-defined electron beam, and this 
again saves considerable weight in the construction of 
the accelerator). Hence, about one-half of the electrons 
remain trapped.* Since the concentration of particles in 
the natural radiation belt is of the order of 10°‘ per 
em’, we can double the radiation intensity after about 
100 seconds of operation. 


After building up the intensity in this manner, it will 
decay away slowly to its original value after the ac- 
celerator is turned off. The trapping volume will change 
from a napkin ring into a ring of more nearly circular 
cross-section (Fig. 5). The main point of the experiment 
now is to study the rate and manner in which the 
radiation intensity falls off after the injection.* 


ype ITE 
= ~ ~ 
\ 
A 
HARD 
Yj RADIATION BELT 


1 10 30 60 100 
ENERGY (MEV) 


Fia. 5. Differential energy spectrum of the radiation belt 
flux. The differential flux is shown as a function of the proton 
velocity 8. Note that the differential flux actually increases 
with 6 until the supply spectrum cuts off, or until the trapping 
properties of the magnetic field deteriorate (indicated by dash- 
line) The corresponding kinetic energies and ranges are shown 
giving the thickness of lead required to stop a proton of 
given energy. At very low energies we find the soft radiation 
belt, its contribution depending on latitude, altitude, and 
probably time. Its uncertain variation is indicated by cross- 
hatching. In any case its energy spectrum will fall off very 
steeply. 


After the electrons are released, we would then check 
on the artificial radiation belt, perhaps once a day, by 
means of simple sounding rockets which traverse the 
region of the belt and determine the profile of electron 
intensity with altitude. From a comparison with the 
theoretically predicted variation one may then expect 
to study the trapping properties of the earth’s magnetic 
field, therefore the mean lifetime of the electrons and 
thereby the density of the atmosphere at various alti- 
tude levels. t 


* The details of the calculation are given in the mathe- 
matical appendix. 

+ We have confined our discussion entirely to an experi- 
ment in which the electrons are released and observed at or 
near the equatorial plane. Quite independently, N. C. Christo- 
filos has considered the injection of electrons for very much 
the same purposes but has elaborated on a somewhat different 
technical approach (29). Basically, however, our ideas are very 
similar. 


The Journal of the Astronautical Sciences 7 


Conclusions 


Modification of the radiation belt by artificial in- 
jection of particles promises to become a valuable tool 
for investigations of upper atmosphere densities and of 
the trapping properties of the earth’s magnetic field. 
This in turn will lead to a better understanding of the 
natural radiation belt and of its connection with mag- 
netic storms and aurora. It should be noted that such 
an experiment of artificial injection of electrons can be 
carried out immediately with available rockets and 
fairly direct modifications of existing instrumentations. 
Once such experiments are carried out, many other 
useful applications may suggest themselves. 


Mathematical Appendix 


We consider injection of 2 Mev electrons near 1000 
km altitude at the equator. We derive the magnetic 
field B = M/r’ and the radius of curvature p as 


PE 
eB 


Numerically B = 0.2 gauss and p = 3.3 X 10° cm. It is 
seen that p < B/grad B = r/3 ~ 2.5 X 10° cm so that 
the condition of the constancy of the magnetic moment 
is probably well satisfied. 

If the line of force intersects a given altitude level at 
a latitude \, then the range of pitch angles of particles 
which stay trapped above this level is 7/2 = ao 2 a, 
where a, = arc sin [cos’\(1 + 3 sin’) ““]. If we set 
this lowest level one scale height Ho below r, then with 
r, = 7.37 X 10°cm and Hy ~ 10’ cm, and using the 
polar equation for a line of force, we find cos’ \ = 1 — 
Ho/r, and >} = 6./7°. This gives for a, = 64°. (The 
fraction of the solid angle filled with trapped particle 
intensity is approximately cosa, = 0.44.) Thus the 
latitude (North-South) dimension will be approxi- 
mately 2rk = 2 X 7.37 X 10° X .116 = 1.71 X 10° cm. 
The radial dimension will be at least ~4p = 1.32 & 10° 
cm, but can be controlled by controlling the time of 
release of the electrons. The azimuthal dimension will 
be 2rr = 4.63 X 10° cm. Thus the volume will be 
~10" cm’ to begin with if we inject over a 10 km 
altitude interval near the equator. 

To calculate the trapping time we set up a diffusion 
equation for the trapped particles: 


an 


AF div 7 + q 


One-dimensional Treatment 


The concentration n, diffusion current 7 and source 
function q are all functions of pitch angle a, energy H 
and time ¢. But in the first instance we will consider F 
fixed and investigate only the diffusion in a. We define 
a pitch angle diffusion coefficient 


p = (Aa! 
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Then we can express 
pe 
Ue 0a 


a an an 
ivi = — —}= D—. | 
LEA abe (v a oar | 
q is assumed to be zero after t = 0, so that finally we 


have dn/at = Da’n/da’. Setting n(a, t) = N (a)-T(t) 
and applying well-known procedures we obtain 


and 


nee ih c(k) eit? ak 


as the most general solution. We define a as 1/2 — 6, 
and allow it to be negative and positive. We now use 
the initial condition: at t = 0, = no(a@) which can be 
expressed as 


role) = / c(k)e** dk 
so that 
if 34 1\ ike’ , 
CUk a= x | no(a’ Je da’. 
20 =) 
Substituting these expressions for c(k) and integrating 
with respect to k we get 
1 2] 


—(a—a’')2/4Dt 
/ no(a’ Je (a—a’)2/ dat 
ioc} 


na, t) = (4xD)2 ¥ 


We take m(a) to be constant over the limits 
—a, <a < a, and zero outside these limits. Then the 
solution can be written 


ta: a&— a ata 

cae {6 Sans | mas | aa | 
es: ie —a |} 
~ © #| appr | J 


where ¢$(x) is the error function, defined as 


9 x 
ped sage 
<a [ (7 dt 


It is useful to examine the behavior of n(0, ¢), i.e., 
the decrease with time in the concentration of particles 
having a pitch angle ay of 90° (and therefore confined to 
the equatorial plane). 


n(0, t) = Coe kes 


When t = 0, [a] = x =, and ¢la,| = 1. At first 
¢ drops rapidly and after a time ty. reaches the value 
0.5. At that point « ~ 0.5 also, so that 


nla, t) 


1 Osean | ee te 
Hb 9) | apn | oe liye D 


For ¢ ~ 0.1, x = 0.1 also; ..t.1 = 100a//4D = 25h 
These results for ti. and to. are fairly representative. 
In greater detail we can consider the variation for other 


values of a, essentially the function 


#| Goon | 


A two-dimensional plot is given in Figure 6. It will be 
seen that the function ¢ decreases very rapidly for large 
values of a (close to a). This expresses the fact that 
after a long time particles will be confined closely to the 
equatorial plane (see Fig. 5). 

_ If we define 


N(t) - | n(a, t) cos ada, 
0 


then N (¢) is given by the area under the various curves 
of Figure 7, weighted by a cosine function. The half 
time for V(t) is seen to be somewhat less than t/2. , but 
fonly by a factor at most. 

- To obtain a rough numerical estimate we use the 
value of D calculated earlier, i.e., 


D=k “ (for relativistic electrons) 

‘where k is 10-20 depending on the nature of the at- 
mosphere. Using y ~ 4, and d ~ 10’ g/cm’* we get 
D~3 X 10‘ rad’/sec. With a, = 64°, the particles 
must diffuse in pitch angle by about 26° or 0.5 rad. 
Therefore t1/2 ~ 10° see or 2 weeks. 


INITIAL VOLUME 
("NAPKIN RING") 


(2) 


(3) 


(b) 


Fic. 6. The initial region into which the electrons diffuse 
is shaped like a napkin ring. It gradually approaches a circular 
cross section as electrons with small pitch angles are lost and 
only electrons with a ~ 90° remain. The exact dependence 
with time can be determined from the solution of a diffusion 

equation. See also Figure 7. Some radial diffusion of the 
electrons will also take place due to small angle scattering. 


$(@t) 


Fig. 7. Solution of the one-dimensional diffusion equation 
indicating the depletion of electrons as a function of their 
pitch angle. 
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Encke’s Method and Variation of Param- 
eters as Applied to Re-entry Trajectories” ' 


Robert M. L. Baker, Jr.t 


Abstract 


The Encke and variation-of-parameters method of pertur- 
bation calculation are developed for the analysis of re-entry 
trajectories. The new techniques have the advantages of both 
higher accuracy and greater computational speed over the 
commonly employed re-entry calculations that involve 
integrations of all the forces acting on the vehicle (known to 
astronomers as Cowell’s method). Applications to both re- 
entry from satellite and interplanetary orbits, which include 
lift as well as drag, will be mentioned. 


Introduction 


The principal deterrent to increased precision in 
trajectory analysis by numerical integration is the 
accumulation of end-figure error. If a problem formula- 
tion is employed that requires the integration of the 
total forces acting on a vehicle (a technique known to 
astronomers as Cowell’s method) large end-figure errors 
may occur, and an additional penalty in calculation 
time is also paid. Even though the availability of large 
scale computing machinery tends to minimize this latter 
factor, the ability to control end-figure error becomes 
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Philosophy at the University of California, Los Angeles. 
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more in doubt as the number of numerical integration 
steps increases. 

When perturbation methods are employed, one inte- 
grates only the departures from an analytic reference 
orbit, and the number of steps required for a given 
accuracy is drastically reduced, with a corresponding 
reduction in both end-figure error and computation 
time. 

Two classes of perturbation methods are commonly 
employed for orbit determination by investigators 
familiar with celestial mechanics. The first of these, 
called ‘‘Encke’s” method, utilizes a reference orbit, and 
deviations from this reference orbit resulting from 
perturbative forces not included in it are accumulated 
during the computational process. If the reference orbit 
selected is a conic section, for example, only the devia- 
tions of the force field from a central field need be 
integrated, and these lead to perturbative displacements 
from the reference conic. 

An alternate approach involves the continuous cor- 
rection (“rectification”) of a reference orbit such that, 
at every instant, the reference orbit exactly conveys 
the position and velocity of the object. This technique 
is designated as the “variation of parameters” method. 
Suitable parameters for the 2-body Keplerian ellipse, 
which can be corrected continuously through the tra- 
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Fig. 1. Perturbation Techniques 


jectory, include eccentricity, semi-major-axis, and lon- 
gitude of perigee. 

The reference orbit that is to be perturbed is usually 
chosen to be the orbit that best depicts the gross mo- 
tion of an object. All forces not taken into account in 
this gross motion are then categorized as perturbative 
forces. It is usually advisable to utilize as the reference 
orbit, an orbit that can be integrated analytically in 
closed form at the outset without requiring a step-by- 
step numerical procedure. These considerations im- 
mediately led astronomers to the utilization of the 
‘two-body Keplerian orbit as the reference orbit for the 
study of planetary perturbations. All forces exerted by 
drag, other planets, etc. are then incorporated as 
perturbative forces that tend to change the orbital 
parameters. 

In the case of re-entry orbits, the drag forces pre- 
dominate, and it is therefore advisable to discard the 
two-body orbit and substitute in its stead a purely 
gravity-free, drag orbit. Because all the drag forces are 
directed opposite to the velocity vector, the orbit will 
be rectilinear. Fig. 1 indicates the difference between 
the drag-perturbed and gravity-perturbed orbits. 


Encke’s Method 


If we denote & and 7 as the x and y directed depar- 
tures for the reference orbit, P, and P, as the x and y 
directed perturbative forces, and the zero subscript 
as the value of the quantity on the reference orbit, 
then Encke’s method (in two dimensions) can be 
formulated (cf. (1)) as 


r= [@-%) + Plat (1a) 


n= | [Gy ) + Plat (1b) 


where 4) and % represent the « and y components of 
acceleration that the re-entry body would have at 
time ¢, if it were moving on the reference orbit. Con- 
sider the coordinate system shown in Fig. 2. The rather 
unconventional labeling of the coordinate axes is em- 
ployed in order to preserve the direction of the x-axis 
towards perigee. 


GRAVITY-FREE TRAJECTORY 


FIXED ARBITRARY DIRECTION 
IN SPACE (i.e., SPECIFIES 
AN INERTIAL REFERENCE FRAME) 


Fig. 2. Coordinate System 


Letting Ap represent the true drag acceleration, we 
find that 


t= Ap, i= A (2a, b) 
a0, Go = (1/Ho) (dijo/dyo), (3a, b) 
P, = —2/r + #4. + ¢,', and (4a) 
Dek Gag OO UA ete ee (4) 


where the A subscripts stand for perturbation that 
result from the asphericity of the Earth and the L 
subscripts stand for lift perturbations. 

Since, in the reference orbit, all of the motion of the 
vehicle is directed along the y-axis, the basic drag 
equation of motion in this orbit can be given by 


i = -| Ao po Vio 


6 = ay) 2 a oA 
eT | TYMAYo ooyHayy (5) 


where Cp, is the reference value of the drag coefficient, 
Ao is the initial projected frontal area of the vehicle, 
mo is the initial mass of the vehicle, go is the acceleration 
of gravity at unit distance (979.81979 cm/sec’), po is 
the sea-level atmospheric density, oc = p/p, 
vy = Cp/Co,, » = m/m, a = A/Ay, and ¥ is the 
inertial speed of the vehicle with respect to the (fixed) 
atmosphere. Here p, Cp, A, and m respectively repre- 
sent the true value of the atmospheric density, drag 
coefficient, area, and mass of the vehicle, while Vo 
is the surface circular-satellite speed (7.905 km/sec). 
Actually, y can be expected to be a function primarily 
of the three variables o, y and 7, , the surface tempera- 
ture of the vehicle. These variables will determine the 
transitional and continuum variation of the drag co- 
efficient (2). In certain re-entry problems in which 
the transitional variation and compressibility variation 
of Cp (and, hence, of y) occur at different points along 
the trajectory, it is possible to separate y into two 
factors: 


Y= yo) xu) (6) 


The function y(c) specifies the variation of the drag 
coefficient during the transition from free-molecule to 
continuum flow, and is strongly a function of the 
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local molecular mean-free path; consequently, of the 
density ratio o. The function y(y) specifies the com- 
pressibility variation of Cp with the Mach number; 
hence, is properly a function of y. If it is assumed that 
ua remains constant and equal to unity during re-entry* 
on the reference orbit and that the gross-variation of 
cy(c) can be obtained from a table or analytically 
by means of a power series, while y(y) can be approxi- 
mately represented by 


y(y) =A + B/y, (7) 
then 
Go = dyjo/dt = —Dy [A + B/yoly(c)oyo (8) 


The reference orbit can be computed in advance either 
numerically (or analytically through employment of 
Eq. (17)) and the value of Eq. (8) is predetermined 
for every point. 

In the true orbit Ap must be computed by employ- 
ing the exact variations of o, y, » and a according to 


Ap = — Di oyna’, (9) 


where v is the velocity of the re-entry body with respect 
to the atmosphere and, hence, must include the effects 
of cross-winds and the rotation of the Earth’s at- 
mosphere. t 

Thus, equations la and 1b can be integrated and 
Encke’s method carried out. 


Variation of Parameters 


In the case of the variation-of-parameters method 
discussed in (5), it is necessary to consider the recti- 
linear orbit only with respect to a spherically sym- 
metrical Earth; therefore, the one angle, #, and one 
distance, x, are the only parameters needed to specify 
the orbit. (Of course, considerations of the asphericity 
of the Earth and rotation of the atmosphere can be 
included also in the perturbation analysis.) Another 
quantity (e.g., y or v) is required to give the position 
of the object on this orbit. The orbital parameters, 
perigee distance, and the longitude of perigee, a, vary 
only because of gravitational perturbative forces. Lift 
forces, normal to the trajectory, also can be easily 
included in the new perturbation scheme. 

Let the grave symbol (‘) denote the perturbative 
variation of a quantity with respect to time in which 
the variation results from gravity. Thus, the radial 
and tangential perturbative accelerations are: 

re ria () (10) 
where r is the radial distance and v is the true anomaly, 
i.e., the angle between r and the z-axis. From (5), we 


* The variation of ua has been included, however, in (3). 

} For example, in equatorial XY, Y, Z coordinates, 
vx =X + YO6,vy = Y — XO, vz = Z, and @ = 27 (1.002,737,- 
803) Harth radii/mean solar day. 
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determine the equations that are to be integrated a 


follows: 
a= fo ‘dt, (11) 
t= / x ‘dt, (12) 
y= it fy + y'| dt, and (13) 
y = i ly + y'] dt. (14) 


If cy(c) and y(y), respectively, can be represented by 


oy(a) = Co + Cor? + Cur* + aes (15a) 
de Vg 45 Yoy" 25 Yu" Spe (15b) 
and yy) = A BF y ees (16) 

then we can establish the equation for | y dt as 
[aa ty + BA) = 
af , 17 

re In ee Tiles y ] 
: = a — B/A. 
exp | ADY 2. FT|, 7 


Thus, given y (and the initial conditions y; and y;), 


/ ydt can be explicitly determined in closed form. 


As derived in (5), 


a’ = ay/yr’, (18) 
é = —2/r'*, (19) 
G-= syle (20) 


We add to the analytical integral, / ydt, the numeri- 


cally integrated value of / y‘dt. It should be noted that 


y. can include also differences between the gross varia- 
tion of, e.g., y(y) and its true variation, ie., y(y) — 
(A + B/y). In this sense, 7‘ will be equal to —y/r’ 


plus other perturbative quantities. The sum of ih ydt 


and iF y‘dt is then the true y. In order to obtain y, we 


must add y‘ to y and then integrate, i.e., 


y =f Wy —a°/gr'l at (21) 


Following this procedure, the variation-of-parameters 


method can be carried out. 


Applications and Conclusions 


The perturbation and computational approaches of 
the foregoing sections are justified only if they lead to 


faster and more accurate calculations. Consequently, 


other calculations were performed in which the con- 


ventional integration of all the forces (Cowell’s method ) 


i 


) 


| 


and the variation-of-parameters scheme were em- 
ployed on the same problem. The problem chosen for 
analysis was the return to the Earth from an inter- 
planetary journey and the Runge-Kutta integrational 
method coupled with the E. Adams correction pro- 
cedure (4) was employed in the actual computations. 
As discussed in (5), the conventional calculations 
(even though their equations were slightly simpler 


and more symmetric) required from 3 to 10 times more 


computer time than did the variation-of-parameters 
computations; furthermore, the many extra integra- 
tional steps involved in the non-perturbation methods 
caused a gradual accumulation of end-figure error. 
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Abstract 


The variational problem of most efficient utilization of the 
fuel is considered for the case of a rocket with separate fuel 
and energy sources and which operates in field-free space. It 
is found that the exhaust velocity should be low (but not zero) 


} initially and should rise, during powered flight, in a manner 
| described in the text. This solution leads to a natural definition 


of propulsion efficiency for such a system which is invariant 
to choice of inertial reference frame. Some properties of the 
variational solution are discussed briefly. 


Introduction 


“Propulsion efficiency” is a term which does not 
appear at present to have a well-defined meaning or to 
be well understood when applied to the topic of rocket 
propulsion. For present-day chemical rockets, the term 
has little significance since maximum performance 
usually results by maximizing the energy content per 
unit mass of the propellant if the density is reasonably 
high and its chemical properties render it compatible 
with the operational requirements of the system. ‘This 
situation results from the fact that the energy available 
resides within the propellant and has an upper limit set 
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Optimum Exhaust Velocity Programming 
and Propulsion Efficiency 


Robert Fox} 


by the physics of molecules. Also, there is an upper 
limit to the mass ratio of such rockets since the propel- 
lant is carried in a container. 

Propulsion systems applicable to travel between the 
planets may not have these limitations. An example is 
an ion propulsion system utilizing nuclear power and 
cesium propellant. Here the energy source is inde- 
pendent of the propellant. On the other hand, even 
though nuclear energy is being used the energy source 
is not infinite and one must consider the problem of 
using the device in the most efficient manner possible. 
Thus we are led to the consideration of ‘‘propulsion 
efficiency”’. 

We shall find that the problem is a dual one. There is 
the fundamental problem of distributing the available 
energy over the propellant mass (variation of exhaust 
velocity) in such a way as to maximize the total 
velocity change. We shall assume that the exhaust 
velocity can be varied without affecting the propulsion 
system mass. In addition, if we are to talk about ‘‘propul- 
sion efficiency’? we must find a unique and meaningful 
definition for this quantity. 


Optimization of the Exhaust Velocity 


Let us consider first the problem of variation of the 
exhaust velocity in the optimum manner for the ideal 
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case of a rocket operating in a space with no gravita- 
tional field. As we shall see, the “obvious” solution, 
letting the velocity v of the exhaust gases with respect 
to the rocket be equal to the instantaneous rocket 
velocity V, is not, in general, the correct one. 

Let the initial mass M, of the rocket under consider- 
ation be composed of a fuel propellant mass M; and a 
dead mass M, which includes the propulsion machinery, 
power plant, control equipment, and the like. Let the 
initial velocity be V; in some inertial frame of reference. 
Let E be the total available energy which can be given 
to the exhaust gases. This is less than the thermal 
energy produced by the power plant since most of the 
thermal energy ends up as rejected heat energy in the 
power plant and propulsion system. We want to pro- 
gram the exhaust velocity v (measured relative to the 
rocket) according to the mass of propellant m which 
has already been consumed. Then, the momentum 
equation is 


(M, — m) dV = v(m) dm 
The final rocket velocity V2 is 


M,=M.+ M;. (1) 


My 
Vi = vit fav = one |b aL (2) 
0 1 = 710 
v(m) is still an arbitrary function and V2 is a function 
of this function. We want to find the maximum value 


of V2 subject to the auxiliary condition 
ny 
pap = ‘| v(m) dm = constant. (3) 
0 


This expresses the limitation on available energy. We 
have implied that the efficiency of the power plant and 
propulsion equipment for converting thermal energy to 
exhaust gas kinetic energy is independent of v. The 
solution of the variational problem posed by equations 
(1), (2), is obtained by solving the Euler equation’ 


ad v(m) 
Ov M, — m 


5, in) 10; (4) 
Ov 


\ is an undetermined parameter which may be evalu- 
ated by using (8). From (4) we have 


= = 1 [4 
v(m) Dr (My ss 7m) . (5) 
Substituting (5) into (3) we find for }, 
A= -— 4/ MT ee 
8EM,(M, — M,)’ 
so that 
__ 1 SENG =H) 
cin ae / Tea 
(6) 


ech 1 2H Moy 
M, Seif) ji 1 


This is the velocity program which maximizes the 
final velocity under the conditions outlined. If we de- 
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note by v the constant exhaust velocity obtained by 
distributing the energy / uniformly over the propellant 


mass 
as Vee 
CSN MY 


then (6) may be expressed 


o( ) Vp ¢ i= . : t 
where the mass ratio iM is 
= : 8 
M ] é . ( ) 


Note that the initial and final values of v(m) are 


v(M;) = Sud. (9) 


(0) =, 


Also note from the very nature of equations (2), and 
(3) that the result is invariant to a transformation 
from one inertial frame to another since only V2, V1 are 
affected by such a transformation and the change in 
rocket velocity V2 — Vi is invariant. V does not appear 
inside the integrals. Substituting (6) into (2) we have 
for the velocity change 


Sis 2EM; — _ , /2H&=1)' Goy4 
Vs Ved ae M2 oe 


since M, = M, — M,. This may be rewritten 


3 (Va— Vi) = #= 48 (11) 
if the exhaust velocity v is programmed according to 
(6). Any other exhaust velocity program would result 
either in reduced velocity change or in different EL, u 
(propellant mass) requirements. For example, if we 
choose instead v(m) = 3 = constant, then 


Ve — Vi = oIng, 


and the relation between M», (V2 — V1), », E is 


1 Spe ee ele 
V.—Vi) = ifn 2 
5) (V2 1) ae (123 


In this case, (V2 — Vj) is less than that given by (11) 
except in the trivial case » = 1. 

Now let us investigate the connection between the 
velocity program (6) and the one obtained by letting 
the exhaust velocity equal the rocket velocity. In this 
case the exhaust gas velocity as measured in the same 
inertial frame as the rocket velocity is zero, since we 
are assuming that the exhaust gases are directed toward 
the rear. Intuition tells us that this ought to be the 
best solution, insofar as utilization of the energy is con- 
cerned. The rocket equation now becomes 


dV dV da 
v Vo © Ai in (13) 


The solution is 
na 
Le ML. (14) 


Note the singularity for V; = 0. 
Let us also calculate the energy consumed. We have: 


My ' 
p=3/ v(m) dm. 
0 


i} Now, 
Vim) M, 
} v(m) (m) and V, eee from (14), 
| Therefore, 
ope TL dm eee oe 
) t Gs 0 (M, <, m)? 2 
3 (15) 
i a8 M.V, R= 1 
{ 2 m 
‘Wor, since V. = pVi, 
t 
2 
I E = (V2 aes Vi) lid M, : (16) 
| 2, wp—l 


We see that this equation is identical with (10) from 
‘the variational solution. There is, however, a distinc- 
| tion. Here we have the auxiliary condition V2/V; = u 
i}from (14), so that for a given value of the initial 
| velocity V; and the mass ratio yw, both # and V2 are 
|} determined. In the variational case H was arbitrary. 
‘| Thus, the velocity program v = V, with uw = V2/V; as 
an auxiliary condition, is a particular solution of the 
variation problem but is not the general solution for 
arbitrary #. 

A further difficulty with the program v = V is that 
it depends upon the choice of inertial frame in which V 
is measured. Observers in two different inertial frames 
a and b could not simultaneously apply such a program 
} to a given rocket. Using these two observers, we can 
| now elucidate further the connection between the 
} variational solution and the program v = V. 

Consider the variational problem for the case V;" = 0 
as measured by an observer in an inertial frame a. 
Let the mass ratio be w and the available energy LH. 
Now transform the observer (not the rocket) to an 
inertial frame b so that the observer then sees 


Ve 
View” 
The relative velocity of the two inertial frames is 
17 a) a b Veo" 
Ue tose 4A Oreste (17) 


since V," — Vi" = V2’ — Vi’. 
In frame a the initial value of the exhaust velocity is, 
from (6), 


“(0) = 1 2EM,(M, — M;) 
aoe ate a M; 


Substituting for Z from (10), 


a M. i M fi a V2" 
= ——— 
v"(0) My, Ve al 
The value of the exhaust velocity as measured by the 
observer in frame b is then 


v(0) = v°(0) — V = 0. 


Furthermore, by substituting (6) into (1) it can 
readily be shown that the value of the exhaust velocity 
seen by an observer in frame a (not an observer on the 
rocket) is constant during the entire burning period. 
Thus, the variational solution for optimum use of the 
energy and starting from rest in frame a is physically 
identical to the solution » = V in the other inertial 
frame b which is related to a through (17). 


Propulsion Efficiency 


We are now in a position to discuss the concept of 
propulsion efficiency as applied to rockets. To begin, 
we might ask what properties such a quantity should 
have. It has already been shown how one can obtain the 
maximum velocity change in a rocket for a given mass 
ratio and available energy. Clearly, since the velocity 
change one can obtain in a rocket is the prime index of 
its performance (for missions in field-free space), this 
solution should correspond to maximum propulsion 
efficiency. Further, since the solution we obtained is in- 
variant to a transformation from one inertial frame to 
another, the propulsion efficiency must also have this 
property. Third, the propulsion efficiency should never 
be negative or exceed unity. To summarize, using the 
symbol 7 to designate propulsion efficiency: 

1. Maximum 7 corresponds to the variational so- 
lution (6). 

2. 7 is invariant to translation of the observer from 
one inertial frame to another. 

oe OuSen seb 

From equation (10) we can immediately find a 
quantity which satisfies all these requirements, namely, 


—_ M.(V2 — Vi)” 
M oF 


(18) 


For the optimum exhaust velocity program, we have 
from (10) 


ae ee (19) 


Clearly, with this definition 7 S 1, since this exhaust 
velocity program maximized (V2 — V,), and since M, 
and # are fixed quantities. Further, 7 2 O for any 
exhaust velocity program since all the quantities in- 
volved are positive definite. The invariance is obvious. 
Also, the problem of finding an exhaust velocity pro- 
gram that maximizes 7 is identical with the variational 
problem stated earlier. 

It is interesting to consider three other possible defi- 
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nitions for 7, which shall be designated m , 72, and 73 , 
and investigate the conditions under which they fail. 


Ey’ — B M, Vi — M2 V2 
1 air aa = 1( 228 E ae (20) 
M.(V2 — Vr) 
ee pees Seer ee PAL 
ne 35 (21) 
ee rate of doing work (22) 


rate of total energy loss’ 


Intuitively, (20) would appear to be more reasonable 
than (18). In fact, if V; = 0, conditions 1, 3 above are 
satisfied. However, condition 2 is not. Further, if 
V, = 0, condition 3 is not always satisfied. For ex- 
ample, consider a rocket with constant exhaust velocity 
v. Then, 


VS Via pins (23) 


LOS Gey 
an <(F)- a 


Substituting (24) into (20) we have 


Let 


m = 0 
for such cases. 

It may be observed that if V; = O in an inertial 
frame b and if F is related to M,, Vi, u by (15) so that 
v = V is the solution of the variational problem in this 
frame, then we find 7, = 1. This is to be expected since 
the exhaust gases have zero kinetic energy in this frame 
of reference. 

The definition 72 does not have the deficiency of be- 
coming negative. However, for the particular case just 
discussed we find that 7m. > 1 since the negative term 
is reduced by the factor 1/u so that 


nm > m = 1. 


Finally, let us consider the definition which may be 
: . 2 
found in certain textbooks’, namely, 73 . 
We have, for an element dm of propellant expended, 


dm 


= os dt i, ; ae Vodm Me 
™B1dm 2, 1dm 2 3(dm* + dmv’) (25) 
=—-v +=—_V 
2 dt Be oli 


{ 
The numerator is just force X velocity X time ele- 
ment. The two terms in the denominator are the energy 
given to the propellant and the diminution of rocket 
kinetic energy due to the reduction in propeliant mass. 
Note that the three terms all have different transforma- 
tion properties. Clearly, 7; < 0 in any reference frame 
in which V < 0. 


Conclusions 


It has been shown that for the case of no gravitational 
field the performance of a rocket can be improved by 
varying the rocket exhaust velocity in the appropriate 
manner during the burning period if the energy source 
and propellant source are distinct from one another. 
Note that the rate of burning, or the thrust, does not 
enter into the discussion. The problem of optimum 
thrust programming for real missions in a gravitional 
field is a quite different and considerably more difficult 
problem than the idealized one considered here. The 
foregoing discussion does indicate, however, that for 
missions involving propulsion devices such as ion 
rockets it may be worth while to consider varying the 
exhaust velocity as well as the thrust. For example, if 
one is going “uphill” against the sun’s gravitational 
field and the propulsion system is strongly power- 
limited, the velocity program (6) helps in two ways. 
First, since the initial exhaust velocity is lower, the 
initial thrust can be higher, resulting in greater initial 
acceleration and reduced flight time. Second, because of 
the lower initial exhaust velocity and the power limita- 
tion, the initial propellant consumption rate is also 
higher so that on the average the propellant is dumped 
at a lower gravitational potential. 
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Format of Technical Papers 


for AAS Meetings 


At the suggestion of various members of the 
Society the Technical Papers Committee Chair- 
man, Dr. Horace Jacobs, has prepared a sample 
sheet of the basic requirements for technical 
papers to be submitted for possible inclusion in 
various AAS meetings. This is reproduced at this 
time for wide distribution. Abstracts of such com- 
munications should be typed single-spaced in a 
text width of 4.5 inches centered below the title 
and author(s) of the paper. The abstract can be 
brief and should not exceed 300 words in length. 
1. Introduction. A technical paper for publication 

in the AAS Proceedings should be typed on a 

good quality bond, 83 x 11 inches. Except for 
the abstract, page width is 63 inches and page 
length is 9 inches. Copy should be typed in 
standard elite, preferably on an electric type- 
writer, and should be suitable for offset. Cor- 
rections should be stripped in rather than 
erased. Text is single-space with double-space 
between headings, equations, items, and para- 
graphs. Major heads are all caps and flush left. 

2. Other Main Headings. The paper can be divided 
into principal sections as appropriate. Headings 
or paragraphs are not numbered. Secondary 
headings are in caps and lower case; they are 
flush left and underlined. 

a. Equations. Equations are as follows: 


(a?mo/ ats?) (a?mo/as2?) = 0 (1) 


Slashes should be used to separate numera- 
tor and denominator if it is conveniently 
possible. Equations are referenced in the 
text as equation (1), (2), ete. 

b. Symbols. Slashes are also preferred to 
separate numerator and denominator in the 
text: m; = mo exp (—v:/c). Where super- 
scripts and subscripts make single-space 
typing impractical, the lines where they 
occur can be typed in space-and-a-half 
or in double-space. 

c. Greek Symbols. If tertiary headings are 
necessary they are flush left, upper and lower 
case, and underlined; text is run-in. Greek 
symbols can be prepared in one of the follow- 
ing ways: 

(1) Typed by varitype 


(2) Typed by a mathematical typewriter 

(3) Hand-drawn. If this method is used the 
following should be noted: Use black or 
red ink, Draw symbols neatly. 
Above is an example of format for enumer- 
ations. They are indented 5 spaces. 

d. Tables. Table headings are centered above 

the table and have the following format: 


Table I 
THE HEADING IN ALL CAPS 


e. Figures. Line drawings suitable for reproduc- 
tion may be inserted by the author at 
appropriate places in the text. If this is done, 
figures should be reduced consistent with 
clarity and space conservation. Captions 
are centered below the figure and have the 
following format: 


Fig. 1. The Caption of the Figure 

References to figures are made in the text 
by the designation Fig. 1, Fig. 2, ete. 
If line drawings are not inserted in the text 
they should be provided as original artwork 
or photostats suitable for reproduction. 
Glossy black and white prints should be 
provided for half-tone photographs. Each 
figure should be numbered and a caption 
provided. 

f. References and Footnotes. References listed 
at the end of the paper are indicated in the 
text by a superscript Arabic number’. 
Footnotes are called out in the text by super- 
script symbols *, **, }, t, etc. Standard 
practices footnotes and references are 
acceptable. 

3. Conclusion. The AAS would appreciate it if 
each author would provide the Society with 
at least 100 preprints of his technical paper 
before the beginning of the session. These 
preprints should be presented in the format 
indicated above. However, illustrations 
might be interspersed in the text or placed 
at the end of the paper. 

4, References. 

1. C. S. Draper, Education in the Astronauti- 
cal Sciences, J. Astronautics, Vol. IV, No. 
2, pp. 29-30, 1957 


